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Abstract Under the standard model extension (SME) framework, Lorentz invariance is tested 
in five binary pulsars: PSR J0737-3039, PSR B1534+12, PSR J1756-2251, PSR B1913+16 
and PSR B2 127+1 1C. By analyzing the advance of periastron, we obtain the constraints on a 
dimensionless combination of SME parameters that is sensitive to timing observations. The 
results imply no evidence for the break of Lorentz invariance at 10~ 10 level, one order of 
magnitude larger than previous estimation. 
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1 INTRODUCTION 

Unification of general relativity (GR) and quantum mechanics is a grand challenge in the fundamental 
physics. Some can didates of a self-consistent quantum theory of gravity emerge from tiny violations of 



Lorentz symmetry ( IKosteleckv 



2005 



Mattinglv 



20051) . To describe observable effects of the violations. 



effective field theories could be a theoretical framework for tests. 

The standard model extension (SME) is one of those effective theories. It includes the Lagrange densities 
for GR and the standard mo del for particle physics and allows possible breaking of Lorentz symmetry 



(IBailey & Kostelecky 



20061). The SME parameters s^ v control the leading signals of Lorentz violation 
in the gravitational experiments in the case of the p ure-gravity sector of the minimal SME. By analyzing 
archival lunar laser ranging data, 



Battat et al 



(2007) constrain these dimensionless parameters at the range 



from 10 -11 to 10~ 6 , which means no evidence for Lorentz violation at the same level. 

However, tighter constraints on s M1/ would be hard to obtain in the solar system because the gravitational 
field is weak there. Thus, for this purpose, binary pulsars provide a good opportunity. Because of their 
stronger gravitational fields, for example the relativistic periastron advance in the double pulsars could 
exceed the corresponding value for Mercury by a factor of ~ 10 5 , these systems are taken as an ideal 
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and clean te st-bed for testing GR, alternative relativistic theories of gravity and modif i ed gravity, such as 



the works by 



and 



Deng 



Bell et al 



(1996) 



Damour & Esposito-Farese 



(1996) 



Kramer et al 



(2006) 



Deng et al 



(2009) 



(1201 II) . 

Motivated by this advantage of binary pulsars, we will try to test Lorentz invariance under the SME 
framework with five binary pulsars: PSR J0737-3039, PSR B1534+12, PSR J1756-2251, PSR B1913+16 
and PSR B2127+11C. In Sec. [2] the orbital dynamics of double pulsars in the SME will be briefed. 
Observational data will be used to constrain the SME parameters in Sec. [3] The conclusions and discus- 
sions will be presented in Sec.|4] 

2 ORBITAL DYNAMICS OF DOUBLE PULSARS IN SME 



When the pure-gravity sector of the minimal SME is considered, it will cause secular evolutions of the 
orbits of double pulsars. Since timing observations of double pulsars could obtain its value very precisely, 
the periastron advance plays a much more important role in c onstraining s^" and, with widely used notations 



in celestial mechanics, it reads dBailev & Kosteleckv 

' duo \ 

/ SME 



2006) 



dt 
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where M = mj + m 2 , 5m = m- 2 — mi (m 2 > mi) and s = 1 — (1 — e 2 
coefficients s. and s.. for Lorentz violation with subscripts P, Q and k are projections of s^ v along the unit 
vectors P, Q and k. The unit vector fc is perpendicular to the or bital plane of the binary puls ars, P points 



5m 2na(e 2 — e 

S QQ) + Z2YU2 S Q 



(1) 



e") 1 / 2 . In this expression, the 



from the focus to the periastron, and Q = k x P. By definitions dBailev & Kosteleckv 



2006), s k = s° j k j , 



s Q = s oj Q j , s kP = .<•"/,•'/'•'. s kQ ee TitfQi, s PP = s ij P l P j and s QQ = s^Q^K However, according 
to Eq. ([TJ , it is easy t o see that the measurement o f u> is sensitive to a combination of s^ v instead of its 
individual components. [ 



Bailey & Kosteleckv (2006) define the combination as 
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(2) 



and crudely estimate its value at the level of 10 -11 . 

Together with the contribution from GR, the total secular periastron advance of a double pulsars system 

is 



W)' <> 



llz 



2/3 



\27ry \M Q 
where T Q ee GMq/c 3 = 4.925490947 /is and 

Ph 
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(3) 



(4) 



The quantity x in Eq. (@]i is the projected semi-major axis, which is usually given by the timing observations, 
while, in some cases, s could be measured directly so that there is no necessity to evaluate it from this 
equation. In this work, Eq. OJ will be taken to find the constraints on s u with timing measurements of 
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Table 1 Timing Parameters of the Double Pulsars. 



PSR 


Pb(d) 


AT (AT©) 


e 


s 




Reference 


J0737-3039 


0.10225156248 


2.58708 


0.0877775 


0.99974 


16.89947(68) 


Kramer et al. (2006) 


B1534+12 


0.420737299122 


2.678428 


0.2736775 


0.975 


1.755789(9) 


Stairs et al. (2002) 


J1756-2251 


0.319633898 


2.574 


0.180567 


0.961 a 


2.585(2) 


Faulkner et al. (2005) 


B1913+16 


0.322997448911 


2.828378 


0.6171334 


0.733650 a 


4.226598(5) 


Weisbere et al. (2010) 


B2127+11C 


0.33528204828 


2.71279 


0.681395 


0.76762 a 


4.4644(1) 


Jacobv et al. (2006) 



"Derived value according to Eq. ©. 



Table 2 Values of s. 





Group I 


Group II 


Predicted sensitivity 








(Bailev & Kosteleckv. 2006) 




(-1.24 ±0.54) x 10" 10 


(-1.42 ±0.75) x 10" 10 
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3 OBSERVATIONAL CONSTRAINTS 



Long-term timing observations can determine the geometrical and physical parame ters of binary pul 



sars very well. A mong them, PSR J0737 -3039 (IKramer et al 



PSR J1756-2251 ( Faulkner et al 



(IJacoby et al 



2006 ), PSRB1534+ 12 dStairs et al 



Weisberg et al 



2002). 



2010) and PSR B2127+11C 



20051) . PSR B1913+16 
20061) are good samples for gravitational tests. Some of their timing parameters are listed 
in the Table Q] In terms of the estimated uncertainties given in parentheses after u>, the data pool is divided 
into two groups: Group I, all the double pulsars are taken; and Group II, including PSR B1913+16, PSR 
B1534+12 and PSR B2127+1 1C, which have the smallest uncertainties. 

By weighted least square method, the parameter s u is estimated (see Table |2}. The estimation made 



by Group I is s 
comparison 



1.24 ± 54) x 10" 10 and Group II gives 



= (-1.42 ± 0.75) x 10~ 10 . For 
Bailey & Kosteleckyl (120061) propose the attainable experimental sensitivity of s w is 10~ u , 



which is 10 times less than the results we obtain. 



4 CONCLUSIONS AND DISCUSSION 



In this work, we test Lorentz violation with five binary pulsars under the framework of standard model ex 
tension. It finds that s u , which is a dimensionless combination of SME parameters, is at the order of 10~ 10 
whether all five systems are taken or top three systems with the smallest estimated unc ertainties of periastron 



advan ces are used. This value, one order of magnitude greater than the estimation by 
(120061) . implies no evidence for the break of Lorentz invariance at 10 -10 level. 



Bailev & Kosteleckv 



Nevertheless, as mentioned bv lBailev & Kosteleckyl ( 120061) . the secu lar evolution of the eccentri city of 



the double pulsars should be included in the analysis. Its contribution is (IBailev & Kosteleckv 

'de\ 1 



2006) 



dt 



where 



r n(l - e 2 ) 1 / V - 2e)s 
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s e is a combination of coefficients in s^" and sensitive to observations. However, there is lacking of timing 
observations on double pulsars running for a long enough time so that rare observations could show the 
secular change of e. Even though a few numbers could be derived from data, the uncertainties of them are 
quite larger than those of periastron advances. Timi ng observations usually co uld set the upper bounds only, 



such as |e| < 1. 9 x 1CT 14 s" 1 for PSR B1913+16 (Baylor & Weisberg . 



PSR B1534+12 dStairs et al 



1989) and lei < 3 x lCT 15 s" 1 for 



20021) . Hence, we suppose that, at least in current stage, the constraints made 
by e might be looser and the resulting upper bound is \s e \ < 3 x 10~ 10 . Although it is consistent with the 
values of s w we obtain, the exact value of s e remains unknown. Therefore, unless timing observations could 
provide much more definitive results about e, the secular changes of eccentricity would not impose a tight 
constraint on or combinations of s 1 *". 

Another issue for future work is to constrain the components of s^ v directly with double pulsars. 
However, the choice of reference frame affects the values of these components so that a certain reference 
frame must be specified first and the projections of s^ v will be along its standard unit basis vectors. For ex- 
ample, for comparing the constraints due to double pulsars and lunar laser ranging, s^ u has to be projected 
along the same triad of vectors. It means the unit vectors P, Q and k (see Sec|2]) have to be decomposed in 
terms of these vectors, which requires the geometrical information of the orbit of the double pulsars, such 
as the orbital elements fl and uj. Unfortunately, timing observations are not sensitive to those two elements. 
This makes the components of s^ v hard to access directly for now and demonstrates the advantages and 
availability of s u . 
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